16 Period doubling route to chaos 



We now study the "routes" or "scenarios" towards chaos. 

We ask: How does the transition from periodic to strange attractor occur? 

The question is analogous to the study of phase transitions: How does a solid 
become a melt; or a liquid become a gas? 

We shall see that, just as in the study of phase transitions, there are universal 
ways in which systems become chaotic. 

There are three universal routes: 

• Period doubling 

• Intermittency 

• Quasiperiodicity 

We shall focus the majority of our attention on period doubling. 
16.1 Instability of a limit cycle 

To analyze how a periodic regime may lose its stability, consider the Poincare 
section: 




The periodic regime is linearly unstable if 

\xi — Xq\ < \X2 — X\\ < . . . 
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or 

\Sxi\ < I&C2I < • • • 

Recall that, to first order, a Poincare map T in the neighborhood of xq is 
described by the Floquet matrix 



In a periodic regime, 

x(t + r) = x(t). 

But the mapping T sends 

xo + 5x — > xq + M&c. 
Thus stability depends on the 2 (possibly complex) eigenvalues X{ of M. 

If I A^l > 1, the fixed point is unstable. 

There are three ways in which |Aj| > 1: 

Im Xi 




1. A = 1 + e, £ real, £ > 0. &c is amplified is in the same direction: 

> > 



X 7 ^2 X J X 4 

This transition is associated with Type 1 intermittency. 

2. A = — (1 + e). 5x is amplified in alternating directions: 

< • > 



x 3 Xj x X 2 

This transition is associated with period doubling. 
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3. A = a ± i/3 = (1 + e)e ±J7 . \5x\ is amplified, &c is rotated: 




This transition is associated with quasiperiodicity. 

In each of these cases, nonlinear effects eventually cause the instability to 
saturate. 

Let's look more closely at the second case, A ~ — 1. 
Just before the transition, A = — (1 — e), £ > 0. 

Assume the Poincare section goes through x = —0. Then an initial pertur- 
bation xq is damped with alternating sign: 

• • • • • — 

Xj X 3 X 2 Xq 

Now vary the control parameter such that A = — 1. The iterations no longer 
converge: 

• • 

x 

x 2 

We see that a new cycle has appeared with period twice that of the original 
cycle through x = 0. 

This is a period doubling bifurcation. 
16.2 Logistic map 

We now focus on the simplest possible system that exhibits period doubling. 
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In essence, we set aside n-dimensional (n > 3) trajectories and focus only on 
the Poincare section and the eigenvector whose eigenvalue crosses (—1). 

Thus we look at discrete intervals T, 2T, 3T, . . . and study the iterates of a 
transformation on an axis. 

We therefore study first return maps 

Xk+l = f{xk) 

and shall argue that these maps are highly relevant to n-dimensional flows. 
For clarity, we adopt a biological interpretation. 

Imagine an island with an insect population that breeds in summer and leaves 
egges that hatch the following summer. 

Let Xj = ratio of actual population in jth summer to some reference popula- 
tion. 

Assume that next summer's population is determined by this summer's pop- 
ulation according to 

OC j -|- ^ — f OC j S OC j . 

The term rxj + i represents natural growth; if r > 1 the population grows 
(exponentially) by a factor r each year. 

The term sx 2 - represents a reduction of natural growth due to crowding and 
competition for resources. 

Now rescale Xj — > (r/s)xj. Then 

. i ' v • . ry ry* 

ilj J _|_ ^ / ikj J I kL' j * 

Set r = 4u: 

Xj + i = 4fiXj(l — Xj). 
This is called the logistic map. 
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16.3 Fixed points and stability 

We seek the long-term dependence of Xj on the control parameter ji. Re- 
markably, we shall see that ji plays a role not unlike that of the Rayleigh 
number in thermal convection. 

So that < Xj < 1, we consider the range 

< ji < 1. 

Recall that we have already discussed the graphical interpretation of such 
maps. Below is a sketch for fi = 0.7: 



f(x) 



A 



1 




f(x Q ) =Xj 







x, 







Xj 1-1/(4h) 



1 



X 



The fixed points solve 




x* = and x* = 1 , 

where the second fixed point exists only for fi > 1/4 




\Afi(l - 2x*)\ < 1. 
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The stability condition for x* = is therefore 

t* < 1/4. 

The non-trivial fixed point, x* = 1 — 1/(4//), is stable for 

l/4<//<3/4. 

The long-time behavior of the insect population :cfor0</i<3/4 then looks 
like 



1 




1/4 3/4 H 



16.4 Period doubling bifurcations 

What happens for /i > 3/4? 

At ii = 3/4, x* = 1 — 1/(4//) is marginally stable. Just beyond this point, 
the period of the asymptotic iterates doubles: 

fix) 

A 

j 




jcJ 1 X 
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Let's examine this transition more closely. First, look | ' 
at both f(x) and f 2 (x) = /(/(#)) just before the 
transition, at \i = 0.7. 

• Since f(x) is symmetric about x = 1/2, so is 

• If x* is a fixed point of f(x), x* is also a fixed 
point of f 2 (x). 

We shall see that period doubling depends on the 
relationship of the slope of f 2 (x*) to the slope of 

/(**)• 

Feigenbaum, Fig. 2. 

The two slopes are related by the chain rule. By definition, 

X\ = f(x ), x 2 = f{xi) =^ x 2 = f 2 (x ). 
Using the chain rule, 

/ 2 W = ^ /(/(*)) L 




= f'Mf'ifM) 



= f , (x )f , (x l ) 



Thus, in general, 



f n '(x ) = f'(x )f , (x l )...f , (x n - 1 ). 
Now, suppose xo = x*, a fixed point of /. Then 

X\ = Xq = x* 



and 



/ 2 V) = /V)/V) = I/V)I 
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For the example of ji < 3/4, 

l/V)l<i => l/V)l<i- 

Moreover, if we start at x$ = 1/2, the extremum of /, then equation (32) 
shows that 

/'(l/2) = / 2, (l/2) = 

=4> x = 1/2 is an extremum of / 2 . 

Equation (32) also shows us that f 2 has an extremum at the xq that iterates, 
under /, to x = 1/2. These inverses of x = 1/2 are indicated on the figure 
for /i = 0.7. 



What happens at the transition, where \i = 3/41 
At // = 3/4, 

/V) = -1 / 2 (^*) = I- 
Therefore f 2 (x*) is tangent to the identity map. 

Feigenbaum, Fig. 3, fi = 0.75. 




Just after the transition, where fi > 3/4, the peaks of f 2 
increase, the minimum decreases, and 

|/V)| > 1 ==> \f 2 '(x*)\ > 1. 
/ 2 develops 2 new fixed points, and X2, such that 

We thus find a cycle of period 2. The cycle is stable because 
\f(xl)\<l and |/ 2 '04)|<1. 

Feigenbaum, Fig. 4, fi = 0.785. 

157 




Importantly, the slopes at the fixed points of f 2 are equal: 

f\x\) = f{xl). 

This results trivially from equation (32), since the period-2 oscillation gives 
f\x\) = f\x\) f'(xl) = f'(xl) f\x\) = f\xl). 

In general, if x{, x\, . . . , x* n is a cycle of period n, such that 

x r +i = f( x r )i r = 1, 2, . . . , n — 1 
and x\ = /(x*) 

then each x* is a fixed point of f n : 

x r = f (x r ), v = 1, 2, . . . , n 

and the slopes f n \x* r ) are all equal: 

= f'(x*i) f'(xl) ■ ■ ■ /'«), r = 1, 2, . . . , n. 

This slope equality is a crucial observation: 

• Just as the sole fixed point x* of f(x) gives rise to 2 stable fixed points 
x\ and x* 2 of f 2 (x) as increases past fi = 3/4, both and £2 gi ye r i se 
to 2 stable fixed points of f 4 (x) = f 2 (f 2 (x)) as ^ increases still further. 

• The period doubling bifurcation derives from the equality of the fixed 
points — because each fixed point goes unstable for the same fi. 

We thus perceive a sequence of bifurcations at increasing values of fi. 

At /i = fi\ = 3/4, there is a transition to a cycle of period 2 . 

Eventually, \i = fii, where the 2 1 -cycle is superstable, i.e., 

f\x\) = f\xl) = 0. 
At fi = fi2, the 2-cycle bifurcates to a 2 2 = 4 cycle, and is superstable at 
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We thus perceive the sequence 
where 



• fi n = value of fi at transition to a cycle of period 2 n . 

• ft n = value of fi where 2 n cycle is superstable. 

Note that one of the superstable fixed points is always at x = 1/2. 




Note that in the case fi = ^2, we consider the fundamental function to be ^2, 
and its doubling to be f 4 = f 2 (f 2 ). 

In general, we are concerned with the functional compositions 

/ 2 " + ' = / 2 "(/ 2 ") 
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Cycles of period 2 n+1 are always born from the instability of the fixed points 
of cycles of period 2 n . 

Period doubling occurs ad infinitum. 



16.5 Scaling and universality 

The period-doubling bifurcations obey a precise scaling law. 
Define 

^oo = value of fi when the iterates become aperiodic 

= 0.892486 . . . (obtained numerically, for the logistic map). 

There is geometric convergence: 



A^oo — 5 n for large n 



That is, each increment in fi from one doubling to the next is reduced in size 
by a factor of 1/5, such that 

5 n = ^ n+1 — — ► 5 for large n. 

The truly amazing result, however, is not the scaling law itself, but that 



5 = 4.669 



is universal, valid for any unimodal map with quadratic maximum. 

"Unimodal" simply means that the map goes up and then down. 

The quadratic nature of the maximum means that in a Taylor expansion of 
f{x) about x max , i.e., 

/(^max H~ = /(^max) ~\~ &f (^max) H (^max) H~ • • • 

the leading order nonlinearity is quadratic, i.e., 

/"0w) 7^ o. 
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(There is also a relatively technical requirement that the Schwartzian derivative of / must be 
negative over the entire interval (Schuster)) 



This is an example of universality: if qualitative properties are present to 
enable periodic doubling, then quantitative properties are predetermined. 

Thus we expect that any system — fluids, populations, oscillators, etc. — whose 
dynamics can be approximated by a unimodal map would undergo period 
doubling bifurcations in the same quantitative manner. 

How may we understand the foundations of this universal behavior? 

Recall that 

• the 2 n -cycle generated by / 2 ™ is superstable at ji = p, n ; 

• superstable fixed points always include x = 1/2; and 

• all fixed points have the same slope. 

Therefore an understanding / 2 ™ near its extremum at x = 1/2 will suffice to 
understand the period-doubling cascade. 

To see how this works, consider f^ipc) and f^ 2 (x) (top of Figures 5 and 7). 

The parabolic curve within the dashed (red) square, for f^ 2 {x) looks just like 
fii^x), after 

• reflection through x = 1/2, y = 1/2; and 

• magnification such that the squares are equal size. 



The superposition of the first 5 such functions 
(/> / 2 ? / 4 ? / 8 ? f W ) rapidly converges to a single func- 
tion. 



Feigenbaum, Figure 8. 
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Thus as n increases, a progressively smaller and smaller region near fs max- 
imum becomes relevant — so only the order of the maximum matters. 



The composition of doubled functions therefore has a "stable fixed point" in the space of functions, 
in the infinite period-doubling limit. 

The scale reduction is based only on the functional composition 

r +1 =/ 2 "(/ 2 ") 

which is the same scale factor for each n (n large). 

This scale factor converges to a constant. What is it? 

The bifurcation diagram looks like 

1 



1/2 




m m [l 2 [l 2 



Define d n = distance from x = 1/2 to nearest value of x that appears in the 
superstable 2 n cycle (for ji = p, n ). 

From one doubling to the next, this separation is reduced by the same scale 
factor: 



d 



—a. 



n+l 



The negative sign arises because the adjacent fixed point is alternately greater 
than and less than x = 1/2. 

We shall see that a is also universal: 

a = 2.502... 
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16.6 Universal limit of iterated rescaled /'s 

How may we describe the rescaling by the factor a? 
For fi = d n is the 2 n ~ 1 iterate of x = 1/2, i.e., 

d„ = /?;-(i/2)-i/2. 

For simplicity, shift the x axis so that x = l/2^x = 0. Then 

The observation that, for n ^> 1, 

~ —a =4> lim (— a) n d n+ \ — r n converges. 

d n+ i n^oo 

Stated differently, 

lim (-a) n /? n (0) must exist. 

Our superposition of successive plots of / 2 ™ suggests that this result may be 
generalized to the whole interval. 

Thus a rescaling of the x-axis describes convergence to the limiting function 



9l (x)= lim(-a) n /C 



71— >00 



n t 2 n 

+i 



X 



{-a) n _ ' 

Here the nth interated function has its argument rescaled by l/(—a) n and 
its value magnified by (— a) n . 

The rescaling of the x-axis shows explicitly that only the behavior of /j^ 
near x = is important. 

Thus g\ should be universal for all fs with quadratic maximum. 

• Figure 5 (top), at fii, is g\ for n = 0. 

• Figure 7 (top), at /i2, when rescaled by a, is g\ for n = 1. 
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• g\ for n large looks like (Feigenbaum, Figure 9) 




The function g\ is the universal limit of interated and rescaled fs. Moreover, 
the location of the elements of the doubled cycles (the circulation squares) is 
itself universal. 



16.7 Doubling operator 



We generalize g\ by introducing a family of functions 

x 



q l = lim (-a) n ff 



2 = 0,1, 



Note that 



Qi-i = \im(-a) n fT 



n— >oo 



X 



{-a) n _ 

= lim ( n)( af 
Set m = n — 1. Then 



-1 f2™" 1+1 

-l+i 



(-a) (-a)"- 1 



(33) 
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and 



9i-i = 



lim (-a)(-a)™C +( i j^, (-*Tf^ 



m r2 r 



X 



-a) {—a) 1 



-agi 



9% 



x 



—a 



We thus define the doubling operator T such that 







' x y 




9% ( 





Taking the limit % — > oo, we also define 

flf(ar) = \m\gi{x) 



l— >oo 



= lim(-«)Vf 



n— >oo 



X 



{-a) n _ 

We therefore conclude that g is a fixed point of T: 

x ^ 



g 



—a 



(34) 



g{x) = Tg{x) = -ag 

g{x) is the limit, as n — > oo, of rescaled / 2 ", evaluated for ji 
Whereas g is a fixed point of T, Tgi, where % is finite, interates away from g 
Thus g is an unstable fixed point of T. 



16.8 Computation of a 

To determine a, first write 

S(O) = -as fe(0)] . 
We must set a scale, and therefore set 

g(0) = 1 ^(1) = 
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There is no general theory that can solve equation (34) for g. 

We can however obtain a unique solution for a by specifying the nature 
(order) of e/'s maximum (at zero) and requiring that g(x) be smooth. 

We thus assume a quadratic maximum, and use the short power law expansion 

g{x) = 1 + bx 2 . 

Then, from equation (34), 

bx 2 \ 



g(x) = 1 + bx = — ag 1 + 



a 2 



—a 



bx 2 \ 2 



1 + 6 1+ , 



2b 2 

= -ail + b) - — x 2 + 0(x A ) 
a 



Equating terms, 

a = -, a = —26 

1 + 6' 

which yields, 



-2 ± V12 

6 = ~ —1.366 (neg root for max at x = 0) 

and therefore 

a ~ 2.73, 

which is within 10% of Feigenbaum's a = 2.5028 . . ., obtained by using terms 
up to x 14 . 



16.9 Linearized doubling operator 

We shall see that 5 determines how quickly we move away from g under 
application of the doubling operator T. 

In essence, we shall calculate the eigenvalue that corresponds to instability 
of an unstable fixed point. 
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Thus our first task will be to linearize the doubling operator T. 5 will then 
turn out to be one of its eigenvalues. 

We seek to predict the scaling law 

— r- — 

H"n /^oo , 

now expressed in terms of fii rather than 
We first expand fp(x) around fp^x): 

ffl(x) ^ ffl^ix) + (fi ~ fioo) Sf(x), 

where the incremental change in function space is given by 

dfp(x) 



Sf(x) = 



Now apply the doubling operator T to fp and linearize with respect to Sf: 



Tfp = -afp 



h 



x 



—a 



-a [fpoo + - Voo) Sf] o 



fi 



x 



—a 



+ (fi - Voo) Sf 



x 



—a 



= Tfp x + (fi - fi^Lf^ Sf + 0(5 f) 
where Lf is the linearized doubling operator defined by 



L f Sf = -a\ f 



f 



x 



-a 



Sf 



x 



-a 



f 



x 



-a 



(35) 



The first term on the RHS derives from an expansion like g[f(x) + 5f(x)] ~ g[f(x)] + g'[f(x)]5f(x). 

A second application of the doubling operator yields 

T(T(fp)) = T'fp^ + (fi- fioo)L Thoo L hoo Sf + 0((Sff). 
Therefore n applications of the doubling operator produce 

T n fp = T n fp x + (fi- fi^) Lm-if^ ■ ■ ■ L Uoc Sf + 0((Sff) . (36) 
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For j2 = poo, we expect convergence to the fixed point g(x): 

~ g(x), n ^> 1. 

Substituting g{x) into equation (36) and assuming, similarly, that Lrf-^ — 
T n f-^x)^g{x) + {fi-li 00 )L n g 5f{x), n»l. (37) 
We simplify by introducing the eigenfunctions <\> v and eigenvalues \ v of L g : 

Lg(j) V = v = 1, 2, . . . 

Write 5f as a weighted sum of <\> v \ 

V 

Thus n applications of the linear operator L g may be written as 

Now assume that only one of X v is greater than one: 

Ai > 1, \ v < 1 for v ^ 1. 

(This conjecture, part of the original theory, was later proven.) 

Thus for large n, Ai dominates the sum, yielding the approximation 

L n g 6f -AJci^i, n»l. 

We can now simplify equation (36): 

r n /^(x) = + (p, - /Zoo) • £ n • a • /i(ar), n > 1 

where 

(5 = Ai, a = ci, and h{x) = 0i. 

Now note that when x = and p, = p, n , 

T n UM = 9(0) + (An - Aoo) • 5 n ■ a ■ h(0). 
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Recall that x = is a fixed point of /|" (due to the x-shift). Therefore 



TV,„(0) = (-«)X(0) = 0. 

Recall also that we have scaled g such that g(0) = 1. We thus obtain the 
Feigenbaum scaling law: 



lim (fj, n - ji^) 5 n = 



-1 



a ■ h(0) 



= constant! 



16.10 Computation of 5 

Recall that 5 is the eigenvalue that corresponds to the eigenfunction h(x). 
Then applying the linearized doubling operator (35) to h(x) yields 



L g h(x) = -a jp' g (^Z^j 



h 



x 



+ h 



-a 



x 



9 



-a 



= 5 • h{x). 



Now approximate h{x) by h(0), the first term in a Taylor expansion about 
x = 0. 

Seting x = 0, we obtain 

-a{g' [g(0)} h(0) + h[g(0)}} = 6 ■ h(0). 
Note that the approximation 

h{x) ~ h(0) => %(0)] = h(l) ~ /i(0). 
Thus /i(0) cancels in each term and, recalling that ^(0) = 1, 

-a [g'(l) + 1]=5. (38) 
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To obtain g'(l), differentiate g(x) twice: 




Substitute x = 0. Note that 

g(0) = and /(0) ^ 
because we have assumed a quadratic maximum at x = 0. Then 

9"(0) = — [</(l)fl"(0)] . 
a 

Therefore 

</(!) = -a. 
Substituting into equation (38), we obtain 

6 = a 2 — a . 

This result derives from the crude approximation /i(0) = h(l). Better approximations yield greater 
accuracy (Feigenbaum, 1979).) 

Recall that we previously estimated a ~ 2.73. Substituting that above, we 
obtain 

£ ~ 4.72, 

which is within 1% of the exact value (!) = 4.669 . . .. 



16.11 Comparison to experiments 

We have established the universality of a and 5: 



170 




reduced by 5 



These quantitative results hold if a qualitative condition — the maximum of / 
must be locally quadratic — holds. 

At first glance this result may appear to pertain only to mathematical maps. 
However we have seen that more complicated systems can also behave as if 
they depend on only a few degrees of freedom. Due to dissipation, one may 
expect that a one-dimensional map is contained, so to speak, within them. 

The first experimental verification of this idea was due to Libchaber, in a 
Rayleigh-Benard system. 

As the Rayleigh number increases beyond its critical value, a single convection 
roll develops an oscillatory wave: 




probe 



Ra=Ra c Ra>Ra c 



A probe of temperature X(t) is then oscillatory with frequency fi and period 

1//1. 

Successive increases of Ra then yield a sequence of period doubling bifurca- 
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tions at Rayleigh numbers 

Rai < Ra2 < Ra3 < 
The experimental results are shown in 



BPV, Figure VIII. 13a and VIII. 13b 



Identifying Ra with the control parameter fi in Feigenbaum's theory, Libcha- 
ber found 

5 ~ 4.4 

which is amazingly close to Feigenbaum's prediction, 5 = 4.669 .... 
Such is the power of scaling and universality! 
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